arXiv:1505.02148v2 [cond-mat.soft] 25 Jun 2015 


Fluids confined in wedges and by edges: 

From cluster integrals to thermodynamic properties referred to different regions 

Ignacio Urrutial* 

^Departamento de Fisica de la Materia Condensada, Centro Atomico Constituyentes, 

CNEA, Av.Gral. Paz 1499, 1650 Pcia. de Buenos Aires, Argentina and 

CONICET 

Recently, new insights in the relation between the geometry of the vessel that confines a fluid and 
its thermodynamic properties were traced through the study of cluster integrals for inhomogeneous 
fluids. In this work I analyze the thermodynamic properties of fluids confined in wedges or by 
edges, emphasizing on the question of the region to which these properties refer. In this context, 
the relations between the line-thermodynamic properties referred to different regions are derived as 
analytic functions of the dihedral angle a, for 0 < a < 27r, which enables a unified approach to both 
edges and wedges. As a simple application of these results, I analyze the properties of the confined 
gas in the low-density regime. Finally, using recent analytic results for the second cluster integral 
of the confined hard sphere fluid, the low density behavior of the line thermodynamic properties 
is analytically studied up to order two in the density for 0 < a < 27r and by adopting different 
reference regions. 


I. INTRODUCTION 

The interest on confined inhomogeneous fluids covers 
a large length of scales of the particles size which starts 
at the simplest one-atom per molecule (e.g. the noble 
gases) and goes up to proteins, polymers (including DNA 
molecules), and large colloids.[1-5] The thermodynamic 
properties of these systems are influenced by the geome¬ 
try of the vessel or substrate that constrains the spatial 
region where the molecules of the system are enabled 
to move. Important efforts are continuously devoted to 
reach a detailed description of the response of fluids to 
some simple geometrical constraints like the confinement 
in pores with slit, cylindrical and spherical shapes, as well 
as the case of fluids in contact with planar and curved 
walls. 

This work focuses on fluids confined by open dihedrons 
built by two planar faces that meet in an edge. Previ¬ 
ous studies were dedicated to analyze the adsorption of 
liquid-vapor coexisting phases on edges and wedges,[6-9] 
and also, to the adsorption on corrugated surfaces. [10- 
12] The main characteristic of the thermodynamics of 
fluids confined by edges and wedges is the existence of 
line tensions. This property also characterizes systems 
with two coexisting phases adsorbed on planar substrates 
(sessile drops) and systems with three coexisting phases 
that meet on a common line.[13, 14] 

One of the particularities of the edge/wedge type of 
confinement is that it produces non-trivial spatial iiiho- 
mogeneities of the fluid. As in the case of fluids conhned 
by curved walls, it happens that different points of view 
in the very beginning of the analysis produce dissimilar 
properties.[15] Thus, it is relevant to establish the basis 
that allow us to compare the thermodynamic properties 
found by adopting these different points of view. 
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In this work, I analyze the statistical mechanics and 
thermodynamic properties of a fluid confined in an 
edge/wedge on the basis of the representation of its grand 
potential in powers of the activity. In Sec. II differ¬ 
ent type of edge/wedge confinements are discussed and 
the thermodynamics of the fluid composed by spheri¬ 
cal particles is revisited. There, I analyze the free en¬ 
ergy and related thermodynamic magnitudes of the con¬ 
fined fluid emphasizing on the explicit choice of the ref¬ 
erence region to which system properties refer. Sec. Ill 
describes the functional dependence of the cluster inte¬ 
grals with the measures of the edge/wedge spatial re¬ 
gion and the consequences that follow on system prop¬ 
erties. There, new relations between bulk- surface- and 
line-thermodynamic properties for different reference re¬ 
gions, are shown. They take the form of transforma¬ 
tion laws and apply to any density. Also, the behavior 
of low density gases is discussed. In Sec. IV, this ap¬ 
proach is utilized to derive analytic expressions for the 
thermodynamic properties (pressure, surface tension, line 
tension, surface- and linear- adsorptions) of the confined 
hard sphere fluid up to order two in density. The con¬ 
sequences of adopting different reference regions on the 
thermodynamics of this system are also discussed in this 
section. The expressions obtained of line-tension and 
line-adsorption, show the dependence with the opening 
dihedral angle. Final remarks are presented in Sec. V. 

II. DETAILED DESCRIPTION OF A FLUID IN 
AN EDGE/WEDGE CONFINEMENT 

Let us consider an open system of particles at constant 
temperature T and chemical potential /i, which is con¬ 
fined by two planar walls that intersect in an edge. For 
simplicity we only refer here to spherical particles. The 
walls exert a hard potential 4> (r) that constrains the po¬ 
sition of the center of each particle to a region A with 
dihedral shape (throughout this work the open-dihedron 
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Figure 1. Fluid confined by a hard-wall dihedron. In the 
region A in light-gray (green) particles are free to move while 
the region in white is forbidden. Note that no matter the 
value of the opening angle both light-gray (green) and white 
regions are straight-edge dihedrons. 



Figure 2. Fluid confined by a hard wall dihedron that induce 
an excluded region. Painted with darker gray is the hard 
wall, lighter gray (green) corresponds to the fluid and the 
excluded region is in white, (a) shows the case a < tt while 
(b) corresponds to the case a > tt. The wall-particle hard 
repulsion distance is a 12 and dashed circles represent particles 
at selected positions near the edge/wedge. The arrows show 
characteristic lengths in cr units. 


confinement is produced by the external hard potential 
being the inter-particles potential arbitrary. In summary, 
Figs, la, lb and 2a correspond to a straight-edge con¬ 
finement while Fig. 2b corresponds to a rounded-edge. 

Before analysing the thermodynamic properties of this 
confined fluid it is necessary to adopt a region B as the 
reference region (RR).[15, 16] Note that B may coincide 
or not with A. I wish to underline that in the study of 
confined fluids is crucial to clearly establish the adopted 
RR which fix the position and shape of its boundary the 
surface of tension. This question is as important as to 
establish the system of reference in the study of a me¬ 
chanical system. Hence, 1 adopt as RR the region B 
which nearly follows the shape of A and has measures 
Mg = (y, A, L) (being V, A and L the volume, sur¬ 
face area and length of the edge of B, respectively). A 
detailed analysis of different prescriptions for B is pre¬ 
sented in Sec.111. The grand potential of the confined 
fluid, relative to B, can be written as 


n = -pv + jA + rL, (1) 

where P is the pressure of the fluid, 7 is the wall/fluid sur¬ 
face tension (or surface free-energy), and T the wall/fluid 
line tension (or line free-energy). The mean number of 
particles in the confined system is 


N = -^ 

dpi ’ 

= pV -\- TaA + Ti^L , 


( 2 ) 

( 3 ) 


geometrical shape is referred to as dihedron) being a the 
inner angle between faces (inner to A). I analyze two dif¬ 
ferent types of edge/wedge confinement. Fig.l shows one 
of the edge/wedge confinement considered. There, the 
edge/wedge available region A is defined by the Boltz¬ 
mann factor exp [—(/(r)/fcT] = 0(|r —C|) where k is 
the Boltzmann’s constant, 0 (x) the Heaviside function 
[0 (x) = 1 if X > 0 and zero otherwise], the dihedral 
region which is the complement of A is C = \ A, 

and ]r —C| is the shortest distance between r and C. 
Note that the faces of A meet on a straight line, thus, 
1 call it a straight-edge dihedron. On the other hand, 
one has the confinement defined by the Boltzmann fac¬ 
tor exp [—(j) (r) /kT] = 0 (]r — C| — being C a solid 
dihedral region and a 12 the minimum distance between 
the center of a particle and the solid dihedron. The lat¬ 
ter case, particularly relevant for colloidal particles and 
macromolecules, is drawn in Fig. 2 where the forbid¬ 
den region between A and C is also indicated. Note that 
for 0 < a < TT the region A is a straight-edge dihedron 
(shown in Fig. 2a). On the contrary, for tt < a < 27r 
(shown in Fig. 2b) A is a rounded-edge dihedron (it has 
a curved end-of-fluid surface). This kind of rounded edge 


where p is the mean number density, F^ is the excess 
adsorption per unit area, and F/, is the excess adsorption 
per unit length. Naturally, this kind of linear decomposi¬ 
tion also applies to other magnitudes, such as the entropy 
S' = — Ip, the energy U = Q+TS+pN, and higher order 
derivatives like — N^ = kT^ which describe 

fluctuations. 

From Eq. (1) is clear that once A is fixed, H be¬ 
comes independent of the adopted RR. On the contrary, 
since measures are relative to S, some of the magnitudes 
(P, 7 , T) depend on B. Naturally, the same argument 
shows that N is independent of the adopted RR although 
some of the magnitudes (p,r^,ri) may depend on B, 
and so on. In summary, even when we have an idea of 
the meaning of the magnitudes (P, 7 , T, p, Fa, Fl) that 
allow us to give a name to each one, they were not ap¬ 
propriately defined yet. Indeed, they do not describe the 
pure properties of the confined fluid (as it may be sug¬ 
gested by the adopted names for these magnitudes), but 
they describe the properties of the fluid with regard to 
a given RR in a sense that will be clarified below. Some 
general aspects of the former discussion follow the anal¬ 
ysis of macroscopic systems with coexisting phases done 
in Ref. ]14]. 
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III. CLUSTER INTEGRALS AND 
THERMODYNAMICS 

Let us consider a system of particles interacting 
through a pair potential iIj (r) with finite range. The 
center of these particles is confined by a hard external 
potential to an edge/wedge region A. It was recently 
shown that the cluster integrals of the system take the 
form[17, 18] 

Ti = i\biV — ilttiA + ilciL . (4) 

Here, the i-th cluster integral r* is linear on the extensive 
measures = {V, A, L) that geometrically character¬ 
ize the region A (its volume V, surface area A, and the 
length of its edge L). We say that are the mea¬ 
sures of the system relative to A. Besides, cluster inte¬ 
gral also depends on the opening dihedral angle between 
faces, a. The volume coefficients hi are the well known 
Mayer’s cluster integrals for homogeneous systems and 
the area coefficients were introduced to describe a fluid 
adsorbed on an infinite wall. [19, 20] hi and with i > \ 
depend on (r) but are independent of a, being ti = biV 
with 5i = 1 and U = Zi, the configuration integral of one 
particle. Eq. (4) was originally derived for the case of 
a straight-edge dihedral region[17] A, i.e. the cases de¬ 
scribed in Figs, la, lb and 2a, and was latter extended 
to the rounded-edge dihedron shown in Fig.2b.]16] 

It is well known that in the low density regime or 
gaseous phase the properties of the confined fluid can 
be rigorously written as power series in the activity 
2 = A“^exp(/3/i) (here /3 = 1/kT is the inverse temper¬ 
ature and A the de Broglie's thermal length).[18] The 
Mayer series of the grand potential for the confined fluid 
is given by 

^ = (5) 

i>i *■ 

being its mean number of particles 

^ = ( 6 ) 

*>i 

By replacing Eq. (4) in Eqs. (5) and ( 6 ), one obtains 

/3L! = - (^ V+(J2 ^ " (E ^ 

i>l i>l i>l 

N = - ibiz'"^ V - ia^z^^ A + ic^z^^ L . (8) 

i>l i>l i>l 

Consistently, a similar transformation applies to other 
thermodynamic magnitudes ]S', U, see Eq. (3)]. Eqs. 
(1) and (7) have a similar dependence with the measures. 
Nevertheless, Eq. (1) is in terms of measures Mg that 
correspond to the choice B as RR, while Eq. (7) is in 


terms of measures regarding to A as RR. Of course, the 
same applies to Eqs. (3) and ( 8 ). Eqs. (1-3) show the 
advantages in deriving as a linear function of the mea¬ 
sures of B. If it is possible, one hnd the 2 ; power series 
for the intensive thermodynamic properties P, 7 , T and 
densities p, L^i, Tg, referred to B. 

Eq. (4) transpires lineal algebra concepts. It can be 
seen as a vector in an abstract space with basis of co¬ 
ordinates {V, A, L) and components (ilbi, —i\ai,i'-Ci) but 
also as the inner product (ilbi, —i\ai,i^-Ci) • M _4 between 
row and column vectors that live in dual spaces. These 
analogies flow to Eqs. (7, 8 ) and will be further investi¬ 
gated in the following Secs. Ill A and IIIB for different 
ranges of a. 

The confinement of the systems drawn in Fig. 1 is 
purely characterized by the region A where the density 
distribution is non-null. This density-based choice of B 
will be labeled with a d subindex (d-RR). For the con¬ 
finement shown in Figs, la and lb the unique simple 
choice for RR is A itself, this prevents to analyze them 
from the point of view of the freedom to choose the RR. 

On the other hand, even when the systems shown in 
Fig. 2 can also be analyzed under the same density- 
based B, other RR could be adopted. To analyze this 
problem and the relation between the thermodynamic 
properties obtained under different choices of B, we study 
the edges/wedges with angles 0 < a < tt and tt < a < 2 tt 
separately. It is interesting to note that when the d-RR 
is adopted the system depicted in Fig. la and the system 
drawn in Fig. 2a are identical, and thus, their properties 
are identical too. 

A. Case I (0 < a < tt) 

For the case of the wedge confinement drawn in Fig. 
2a, I study two different choices for the RR that are the 
most natural to be adopted. Under the density-based 
RR that identifies B with A, the measures are Md = 
(Vd, Ad,i), being the t-th cluster integral 

Ti/i\ = biVd - GiAd + CiiP) L = hd- Md , (9) 

where bd = {bi,—ai,Ci) is the vector of coefficients. The 
second simple choice for B is the empty-region (e-RR), 
i.e. B is taken as A joined with the white region in 
Fig. 2a and the measures are Me = {Ve, A^, L). In 
what follows we will use a (see Fig. 2) as the unit 
length. From geometrical considerations it is possible to 
obtain the linear relation between both sets of measures: 
14 = 14 — ^Ae- 1-7 cot^ L for 0 < a < TT, Ad = Ae— csca L 
for 0 < a < ^, and Ad = Ae — cot^ L for ^ < a < tt. 
Thus, under the e-RR choice 

Ti/i\ = biVe - hjAe -I- Ci(Q;) A = be • Me . (10) 

We introduce the matrix Y that transforms between both 
sets of measures 

Md = Y • Me and Me = • Md . ( 11 ) 
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Its expression follows from the relations above Eq. (10) 


/I -i icotf \ 
y = 0 1 -y{a) , 

Vo 0 1 ; 

(12) 

with y{a) = csca if 0 < a < | and y{a) = cot^ if 
^ < a < TT (note that y{a) is a continuous non-derivable 
function at a = ^). Given that H remains unmodified no 
matter which RR is adopted, one finds the linear relation 

between the unknown coefficients be = [bi,— 

and 

the known bd, through the Y matrix 

/ 

be • Me = bd • y • y-l • Md , 


be = bd • y . 

(13) 

Besides, through Eqs. (1-3) Y also transforms the ther¬ 
modynamic properties 

(-/^,7,r), = (-p,7,r)d-y, 

(14) 

(d,r^,rL)g = (p,rA,rL)j • y, 

(15) 


where the label outside brackets shows the adopted ref¬ 
erence region. These relations should be valid even when 
the series expansion in powers of z does not apply, and 
thus, they should apply to any density. Eqs. (14, 15), 
with Y taken from Eq. (12), are one of the main results 
of the current work. They show the transformation law 
between the intensive properties of the confined system 
when different RRs are adopted. 

Now, we turn our attention to Eqs. (1, 7) and (14). 
They show that = P being P the pressure of 

the bulk fluid at the same T and /r. Furthermore, one 
obtains 

7e = 7d + Pn ■ (16) 

Eq. (16) found here for a wedge confinement is a known 
relation for fluids adsorbed on both planar and curved 
walls. [15, 20] The z power series representation of je and 
7 d shows that they are the surface tension of the fluid 
in contact with an infinite planar wall (each one for a 
different RR). The line tension transforms as 

P a 

%=Td-ldy{a)-—cot-, (17) 

To the best of my knowledge it is the first time that Eq. 
(17), which applies to any density, is derived. Turning to 
Eqs. (3, 8 ) and (15), they imply that Pf. = Pd = P with p 
the number density of the bulk fluid (at the same T and 
p). For the surface adsorption one finds, 

(rA)e = (rA)d-p/ 2 , (18) 

which is known to be an exact relation for planar walls. 
Again, based on the z power series one finds that both 
(r^)g and (ryi)d &re the adsorption of the fluid on an 


infinite planar wall (each one for a different RR). Fur¬ 
thermore, one obtains for the excess linear adsorption 

(hi), = (Ti)^ - (T^)^ yia) + ^ cot| . (19) 

Again, this expression applies to any density and it was 
never published before. 

B. Case II {n < a < 2 tv) 

Now, focusing on the case shown in Fig. 2b, I will ana¬ 
lyze two different choices for the region B, which split in 
three different sets of measures that are the most natural 
to adopt. The first choice is the d-RR, which corresponds 
to identify B with A. For this d-RR one can consider 
two different criteria to dehne the measures depending 
on whether A is taken as the area of the planar part 
(Apd) of the surface dB or as its total area. Thus, us¬ 
ing the first criteria the measures are Mdi = (Vd, Apd, L) 
and 

Ti/i\ = hiVd - aiApd + Ci(a) L = bdi.Mdi , (20) 

with the vector of coefficients bdi = {bi, —ai,Ci). How¬ 
ever, if one adopts the second criteria that assumes A as 
the total area of dB it is obtained 

Ti/i\ = biVd - UiAd + Ci(a) L = bd 2 -Md 2 , (21) 

with Md 2 = (Vd, Ad, L) and Apd = Ad — (a — tt) ^L. The 
relationship between both sets of measures is 

Mdi=T-Md 2 , Md 2 = E-^-Mdi, ( 22 ) 

while the vectors of coefficients relate through 

bd2 = bdi • r , (23) 

with 

/I 0 0 \ 

r= 0 1 -i(a-7r) . (24) 

Vo 0 1 J 

The relations between the equations of state {—P,^,T) 

and also (p, rA,rL), in dl-RR and d2-RR are given by 
Eqs. (14) and (15) with the obvious change of labels 
and with Y taken from Eq. (24). Therefore, one finds 
Pd 2 = Pdi = P (with P the bulk pressure), 7 d 2 = 7 di 
(which are equal to the planar-wall surface tension 7 d 
discussed for the case a < tt) and 

rd 2 = Tdi - 7 d (a - ^) /2 . (25) 

Besides, it is obtained pd 2 = Pdi = P, (rA)d 2 = (rA)di 
(which are equal to the planar-wall adsorption (ryi)^ dis¬ 
cussed for the case a < tt) and 

(rL)d2 = (rL)di - (rA)di (« - ^r) /2. (26) 
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It seems that Eqs. (25) and (26), that apply to any 
density, are novel results. 

The other choice for B is the e-RR, that is, the join of 
region A and the white region in Fig. 2b. In this case, 
the measures are Mo = {Ve,Ae,L) and 

Ti/il = btVe - QiAe + Cj(q;) L = be.Mo , (27) 

with Vd = 14 — ^Ae — i (a — tt) L and Apd = ^e- Eqs. 
(22, 23) describe the transformation between both, the 
measures and the coefficients, they remain valid with the 
change of labels d 2 e and for 


Y = 




1 0 

0 I 



(28) 


Yet non-surprising, following Eqs. (14) and (15) with the 
obvious change of labels and taking Y from Eq. (28) one 
obtains Pe = P and pe = p. Furthermore, both 7 e and 
(r/i)g coincide with the planar-wall magnitudes found 
for the case a < tt, and then Eqs. (16, 18) apply for the 
broad range 0 < a < 27r. Finally, one obtains 

re = rdi+P(a-7r)/8, (29) 


gives surface and linear thermodynamic properties pro¬ 
portional to . Using Eqs. (2, 4) and trivial series ma¬ 
nipulation one obtains the power series for z (p) and the 
series representation of the thermodynamic properties in 
powers of p. Up to order p^ it is obtained: /3P = p — b2p^ 
(i.e. the virial series for the bulk gas [18]) and 

/37 = -r^/2, PT=-Tl/2 (32) 

[with = —2a2p^ and El = 2 c 2 p^]. These notable rela¬ 
tions are not well known. They deal with inhomogeneous 
fluids and link linearly an excess free energy (times f3) 
with the corresponding excess adsorption. It is remark¬ 
able that Eq. (32) does not include coefficients related 
to the interparticle potential. Eq. (32) resembles the 
equation of state of the bulk ideal gas, nevertheless, it 
applies to any edge/wedge confined fluid up to order p^. 
As can be easily verified, the use of e-RR provides more 
complex expressions for the surface and linear thermo¬ 
dynamic properties than Eq. (32). In summary, d-RR is 
appropriate to obtain a simple description for the ther¬ 
modynamics of the confined ideal gas and also of any gas 
at low density, but e-RR is not. 


IV. APPLICATION TO HARD SPHERES 


(rL), = (ri)^,-p(a-7r)/8. (30) 

Eqs. (29) and (30) were not published earlier. 

C. Low density 

In this brief digression the confined ideal gas and the 
low density regime of the confined non-ideal gas are 
analyzed. I first concentrate in the d-RR (d-RR for 
0 < a < TT, dl- and d2-RR for tt < a < 2 tt). To obtain 
the properties of the confined ideal gas one truncates all 
the series in Eqs. (6, 7, 8) at the first order in power of 
z. By adopting d-RR the volume 14 is equal to Zi and 
the first cluster integral is n = 14. Thus, (iP = z, p = z, 

7 = r = 0 , and = El = 0 . (31) 

Therefore, under d-RR the confined ideal gas is thor¬ 
oughly described by fiP = p, i.e. the equation of state 
of the bulk ideal gas. Clearly, if we turn to e-RR the 
edge/wedge confined ideal gas has non-null surface- and 
line- free energies. They can be evaluated using Eq. (31) 
and the transformations discussed in Secs. IIIB and 
Iff A. The conclusion is that in order to obtain the sim¬ 
pler expressions for the thermodynamics of the confined 
ideal gas, the d-RR is better than e-RR. 

For the confined non-ideal gas under d-RR the first 
cluster integral remains unmodified in comparison with 
the ideal gas. The second and higher order could be 
calculated by direct integration. Now, the series given in 
Eqs. (5, 6, 7, 8) are truncated at order two in z, which 


Recently, through adopting the d-RR, the low density 
behavior of the hard sphere (HS) confined fluid in an 
edge/wedge cavity was studied using an analytic expres¬ 
sion of C 2 (a).[I 6 ] In this section we compare those prop¬ 
erties with that found by adopting the e-RR. With this 
purpose the natural units for the HS system will be used 
(which is equivalent to set the particles diameter a as 
the unit length). In Ref. [16] was obtained the following 
exact expression 

C 2 (a) = - [1 -I- (tt - a) cot a] (33) 

that applies for 0 < a < tt in the d-RR, while the analytic 
expression for tt < a < 27r in dl-RR is 

C2{a) =-^ {a - tt) + Q , (34) 

with Q = 0.007125 {1 — exp [—2.74 (a — tt)]}. Using the 
known parameters 62 = —27r/3, 02 = — 7 r /8 and Eqs. 
(33, 34) for 02 ( 0 ;) one readily finds the series expansion 
of { 7 ^, 7 , TAaAl] in power of p up to order two, by 
adopting both d-RR and dl-RR. For T and El it gives 
analytically the angular dependence with a (up to order 
two in p). 

Now, we analyze the consequences of choosing a differ¬ 
ent RR on the thermodynamic properties of the confined 
HS system. In particular, novel analytic expressions of 
relevant line-thermodynamic properties for e-RR and d2- 
RR will be derived. Through the use of Eqs. (14, 15), 
the matrices for RR transformation (12, 24, 28), and the 
density power series of Pd, 7 d, 7d, Tdi, (EA)d, (rL)d> 



6 




Figure 3. Surface tension and surface adsorption. Curves 
for adsorption are marked with squares. In continuous lines 
(blue) are plotted the magnitudes in d-RR while dot-dashed 
lines (red) refers to e-RR. 


(rL)di one obtains the series for 7 , T, and in the 
e-RR and d2-RR. For the walFfluid surface tension and 
excess area-adsorption, both up to terms of order 0{p^), 
it is obtained 


/37e 


p 

2 



(rA)e = 


P 

2 



(35) 


Fig. 3 displays the surface tension and surface excess 
area-adsorption by adopting d-RR and e-RR. There one 
can observe the effect of choosing a different RR in the 
properties of the confined HS fluid. For the d-RR one 
finds / 37 d > 0 and (Fyi)^ < 0 , on the other hand for 
the e-RR they yield < 0 and (F. 4 )g > 0. Besides, 
near p > 0 the null slope in /Syu and (F^)^ is apparent 
while Pje and (F^)^ are linear with density. For Te in the 
range 0 < a < tt there are two branches: the first one for 
a < ^ and the second one for a > ^, the corresponding 
expressions are 


/3re = 

PTe = 


p a p^ 
-4 cot- -h — 


/37r \ 

57r 

(t+“J 

cot a -csca 

8 J 


. , . bn a 

1 + (tt — a) cot a -cot— 

^ ^ 8 2 


(36) 


In Fig. 4 it is shown the low density behavior of the 
edge line tension T for d-RR and e-RR. In both cases, T 
is monotonous. For the Md measures T is positive and 
has positive slope. On the contrary, using Me measures 
T is negative and has a negative slope. For both, the 
modulus of the slope decreases with increasing a. The 
obtained expression of (Fl)^ for 0 < a < tt is 




1 , , IStt 

I + (tt — q:) cot a -\ - 3 — y(a) 


( 37 ) 


Figure 4. Line tension of the HS system vs. density for angles 
in the range 0 < a < tt. Continuous lines (blue) plot 7d 
while dot-dashed lines (red) plot Te- The curves correspond 
to a = 7r/4, tt/2, 37r/4 while the arrow points to the direction 
of increasing values of a. 



Figure 5. Line adsorption of the HS system vs. density for 
angles in the range 0 < a < tt, the curves correspond to 
a = Tr/d, it/2, 37r/4. Continuous lines (blue) plot (F l)^ while 
dot-dashed lines (red) plot (Tl)^- See Fig. 4 for details. 


which is non-derivable at a = ^ [see Eq. (12)]. Fig. 5 is 
similar to Fig. 4 but for the linear adsorption F^. When 
the Md measures are considered the linear adsorption is 
monotonous, negative (i.e., there is local desorption) and 
has negative slope which increases with increasing a. On 
the other hand, when Me measures are adopted Fl is not 
monotonous, it is positive (i.e., there is local adsorption) 
and decreases for larger values of a. 

I also present here a similar analysis for the case 
TT < a < 2tt. The results for the line-tension and ex¬ 
cess linear adsorption by adopting e-RR and d2-RR up 
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Figure 6. Line tension of the HS system vs. density for an¬ 
gles in the range tt < a < 27r. The curves correspond to 
a = 37r/2, 77r/4. In long-dashed lines (blue) it is drawn 

Tdi, short-dashed line (green) corresponds to 7d2 and dot- 
dashed line (red) is for 77- The arrows point to the direction 
of increasing values of a. 


Figure 7. Line adsorption of the HS system vs. density for 
angles In the range n < a < 2n. In long-dashed lines (blue) 
it is draw (rL)^]^, short-dashed line (green) corresponds to 
(r_L)d 2 and dot-dashed line (red) is for See Fig. 6 for 

more details. 
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Fig. 6 shows the linear tension of the HS system in the 
case TT < a < 27r. The functions 7di, 7d2 and Te are 
shown for comparison. In all cases, T is monotonous. 
For the Mdi measures T is negative and has negative 
slope which decreases with increasing a. On the other 
hand, for Me and Md 2 measures, T is positive and has 
a positive slope which increases with increasing a. Even, 
T for Md 2 is nearly zero in the adopted scale. From the 
comparison between Fig. 6 and Fig. 4 the inversion of 
the sign of 7^ at a = tt is evident, where the edge/wedge 
disappears. Fig. 7 plots the linear adsorption of the HS 
system for the same three measure sets. In the cases of 
Mdi and Md 2 we observe a monotonous F/,. For the 
Mdi measures F^ is a positive (i.e., there is local ad¬ 
sorption) increasing function and its slope increases with 
increasing a. On the contrary, for Md 2 measures F^ is 
a negative (local desorption) decreasing function and its 
slope decreases with increasing a. Even, F l for Md 2 is 
nearly zero in the adopted scale. For Mo measures F^ is 
not monotonous, is negative (i.e., there is local desorp¬ 
tion), attains its minimum near p ~ 0.17 and decreases 


with increasing values of a. Fig. 7 and Fig. 5 show the 
inversion of the sign of (F i,)^ when both edge and wedge 
disappear at a = tt. 

In the literature, both d-RR and e-RR were used 
to study HS systems confined in cavities with different 
geometries. [21, 22] The behavior found with d-RR (in¬ 
cluding dl-RR and d2-RR) and e-RR, for both line ten¬ 
sion and line free energy, shows that they strongly depend 
on the adopted reference system. Linear thermodynamic 
magnitudes that are less dependent on the adopted refer¬ 
ence region can be found by considering the mean values 
of excess density and excess free energy in a region with 
finite size around the edge. [16] 


V. FINAL REMARKS 

In this work I studied the relations between the ther¬ 
modynamic properties of fluids confined by wedges and 
edges when different RRs are adopted. The analysis was 
based on the activity series expansion of the grand free 
energy for inhomogenous systems, and on the properties 
of its coefficients the cluster integrals. I utilized a simple 
approach that linearly connects the geometric measures: 
volume, surface area and edge length, in the different 
RR that are considered. From that, the law of transfor¬ 
mation of thermodynamic properties between RRs was 
deduced. A similar method was previously used to study 
a system of hard spheres confined by curved walls. The 
method was here refined and can be used to analyze in¬ 
homogenous fluids confined by walls with a variety of 
shapes. Under this non-standard approach I have stud¬ 
ied the dependence of the linear-thermodynamic prop¬ 
erties on the adopted RR along the complete range of 



















dihedral angles 0 < a < 27r. Analytic expressions that 
transform the thermodynamic intensive properties: pres¬ 
sure, surface tension and line tension of the system when 
different RR are adopted were derived for the first time. 
Surface adsorption and line adsorption were also ana¬ 
lyzed in this framework. The relevant results were given 
in Eqs. (17, 19, 25-26, 29-30). Furthermore, the thermo¬ 
dynamic properties of both, the confined ideal gas and 
of the confined real gases at low density (up to order 
p^), were analyzed by adopting different RR. We found 
that the density-based d-RR is advantageous to obtain a 
simpler analytic description of the studied properties. 

Regarding to the confined HS fluid, which is a rele¬ 
vant reference system both for simple and colloidal fluids, 
the dependence of line adsorption and line tension with 
the edge/wedge dihedral angle and density was analyzed. 


We found explicit analytic expressions truncated to order 
two in density that describe these properties for different 
RR and for the complete range 0 < a < 27r. They are 
shown in Eqs. (36, 37, 38, 39). The new results obtained 
for HS complement those recently published. [16] Given 
that these analytic expressions are exact or quasi-exact, 
they constitute well defined references that should enable 
to validate other approximate theories like fundamen¬ 
tal measure density functional approaches to edge/wedge 
confined fluids at low density. 
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